Magnetostatic fields in tubular nanostructures 
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The non-uniform magnetostatic field produced by the equilibrium and non equilibrium magnetic 
states of magnetic nanotubes has been investigated theoretically. We consider magnetic fields pro- 
duced by actual equilibrium states and transverse and vortex domain walls confined within the 
nanostructure. Our calculations allow us to understand the importance of the magnetostatic field 
in nanomagnetism, which is frequently considered as a uniform field. Moreover, our results can be 
used as a basis for future research of others properties, such as the investigation of spin waves when 
domain walls are present, or the motion of a magnetic particle near a magnetic field. 



I. INTRODUCTION 

Structures of nanometric dimensions are strongly de- 
pendent of size and shape mainly due to the fact that 
the characteristic scales of many physical phenomena 
are comparable to the dimensions of current nanostruc- 
tures. Magnetism makes no exception and fundamen- 
tal research of nanomagnets is further fuelled by their 
prospective applications into technological devices^ and 
biomedical applications 3 -^. The implementation of nano- 
magnets into potential devices requires a deep control of 
their physical properties; hence it is very important to 
control the shape and size of the nanoparticles which are 
crucial to determine their magnetic properties. Within 
currently available nanostructures, nanotubes made from 
3d transition metals and their alloys are of particular im- 
portance, since the nanotube topologies provide us of two 
surfaces for modification, allowing the generation of mul- 
tifunctional magnetic nanoparticles. Clearly, a technol- 
ogy capable to modify different surfaces would be highly 
desirable 5 . 

More recently hollow tubular nanostructures have been 
synthesized^— and may be useful for applications in 
biotechnology, because low density magnetic nanotubes 
(MNs) can float in solutions and are more suitable for 
in vivo applications^. In recent years, the magnetism 
of planar and circular nanowires has been intensely in- 
vestigated, while MNs have received less attention, al- 
though MNs do not exhibit magnetic vortex cores or 
Bloch pointsiS, leading to a different behaviour with a 
more controllable reversal process. 

It is well know that changes in the cross section of 
the MNs strongly affect the reversal mechanism and the 
overall magnetic behaviour^ir— . Recent articles have 
shown that MNs present three main equilibrium states: 
a uniform state (U), a mixed state (M) and a vortex state 
(V) . With basis in theoretical models it has been argued 
that for nanotubes with outer radius less than 4L X (where 



L x = \j2Aj [LqMI is the exchange length, M s the satu- 
ration magnetization and A the stiffness constant of the 
magnetic material) the magnetization is almost uniform 
and oriented principally along the cylindrical axis 1 - 4 -. We 
call this configuration uniform state. Also for nanotubes 
with R > 4L X the magnetic equilibrium state correspond 
to the so-called mixed state, which is a mixture of U 
and V states 8 -^ - — . This state presents a uniform mag- 
netization along the middle region of the tube and near 
the lower and upper surfaces the magnetization deviates 
from the uniformity in order to reduce the magnetostatic 
fieldi*. Finally, there is a transition from the M state to 
the vortex state as can be seen from figure 7 ofi 4 -. 

On the other hand, it has been argued that domain wall 
propagation in nanostructures, which can be achieved by 
external fields or spin-currents, is of basic scientific and 
technological interes t 17 ' 18 . More recently, it has been 
shown that the size-dependent reversal process in MNs 
occurs via DW nucleation and subsequent DW propaga- 
tion. The DW microstructure depends on the tube cross 
section and can be a transverse wall (TW) for tube radius 
(R) smaller than a critical radius (R c (/?)), or a vortex wall 
(VW) for R > R c {pj£r22., where /3 = RJR is the ratio 
between the inner (R4) and outer tube radius. The criti- 
cal radius defined above depends on the magnetic mate- 
rial and on the shape factor j3 and then it can be tailored 
manipulating its geometrical parameters. In polycrys- 
talline systems, the critical radius, R c (0), ranges from a 
few nanometers to 20 nm approximately, and therefore, 
since experimentally fabricated MNs satisfy R > R c {f3), 
we can expect that the propagation of a vortex DW be 
the dominant magnetization reversal mechanism for fer- 
romagnetic nanotubes^—. 

Although the magnetostatic field of MNs with uniform 
magnetization has been already investigated 2 ^—, an an- 
alytical model for the magnetostatic field produced by 
a MN with the above mentioned configurations, either 
equilibrium or transient states, has not yet been per- 
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formed. Our intention is to fill this gap by using simple 
magnetization model a 14 ' 19 . Since nanostructures are usu- 
ally polycrystalline, the crystallographic orientations of 
the crystallites are random and, as a consequence, the av- 
erage magnetic anisotropy of the particle is very small. In 
view of that, it will be neglected in our calculations^^!. 
This paper is organized as follows: In Sec. [II] we present 
the magnetization models, including in Sec. Ill Al the 
equilibrium states M and U and the VW, and then in 
Sec. Ill Bl we discuss the case of the TW. Results and 
discussions are summarized in Sec. Mil and conclusions 
in Sec. HV1 



II. 



CALCULATION OF THE MAGNETOSTATIC 
FIELDS 



In this section we present explicit expressions to evalu- 
ate the magnetostatic field for nanotubes with the mag- 
netization states described above. Within the contin- 
uum theory of ferromagnetism^ 8 ., the magnetostatic field 
is given by H = — VZ7, where U is the magnetostatic po- 
tential, which depends on the magnetization itself. The 
basic assumption we address in this paper is that the 
nanotubes we are considering are thin enough to avoid 
any dependence of the magnetization in the radial coor- 
dinate. In what follows we present the magnetization and 
the related magnetostatic field for the equilibrium states 
(U and M) and both DWs structures aforementioned. 



A. Equilibrium states and vortex-like 
configurations 

First, we consider nanotubes with vortex-like configu- 
rations, as the mixed state and the transient state with 
a VW confined within the tube. The magnetization of 
such magnetic configurations can be generally written 
as M(z) = M s m(z) = M s (m^(z)4> + m z (z)z), where 
m z (z) = cosO(z), and m$(z) — sinO(z), with Q(z) 
being the angle between the local magnetization vector 
and the cylindrical z axis. With this rather general form 
for the magnetization we can compute the magnetostatic 
field produced by M(z). Following the calculations ofii 4 -, 
the magnetostatic potential can be written as 

U(p, z) = -± l g(q)Jo(qp) [fs + fv) dq, (1) 



where f s and f v come from the surface and volumetric 
magnetic charges and can be written as 



f. = -m z (0)e-^ +m z (L)e-^- L K 



fv 



dm z {y) e - q \*- y \ dy 
dy 



(2) 



(3) 



In the potential we have also defined: 

g(q) = (R/q) [J^qR) - PJ^qR)] 



where J m (x) are Bessel functions. The magnetostatic 
field associated with the above potential can be written 
as H = —VU(p, z), and then, the cylindrical components 
of the field are Ha, = 0, 



M s f°° 



dz 



fv] dq, 

dfv 

' dz 



(5) 



dq. (6) 



At this point we are in position to calculate the mag- 
netostatic field through the reversal process for MNs, as 
well as to investigate the changes in the magnetostatic 
field when the equilibrium state is properly considered. 
In what follows we will use the equations above to cal- 
culate the magnetostatic field for the M state and the 
vortex DW. 



1. Equilibrium magnetization states 

The arrangement of magnetic moments in a MN with 
small radius can be seen as a uniform state in the mid- 
dle region of the tube, but with vortex-like deviations 
of the magnetization located at the tube ends&i 4 ^— . 
This mixed state has been observed by micromagnetic 
simulation :) 8 ' 15 ' 16 and studied theoretically 14 . It has been 
shown that two equilibrium states can be studied with 
the model for the M state. In other words, the U state 
is a special case of the former, as we will see later. The 
model we use for the M state is given by the function 14 

( 6 {d-z)/d for < z < d 

6(z) = < for d < z < L - A (7) 

[6 L (z-L + X)/X for L - A < z < L 

where d and A are the dimensions of the regions where the 
magnetization deviates from the z-axis. In the model, 6o 
and 6 1, correspond to angles between the magnetization 
and the z-axis evaluated at the tube ends. Note that 
we obtain the uniform magnetization state in the limit 
m z (0) = m z (L) = 1 (#o = 6l = 0) which occurs approxi- 
mately for a radius smaller than ALj^. The parameters 
of the model are such that they minimize the total energy 
for each set of geometrical parameters. For simplicity we 
consider symmetric nanotubes, that is, tubes without dif- 
ference between both tube ends. In those cases the four 
parameters (d, A, 8q, 9l) variational problem is reduced 
to a two parameters variational problem, when d = A 
and 9q = 9l—- Then, f s and /„ become 



/f = cos0o 



,-l\z-t\ 
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(8) 



}dy. (9) 



With these quantities in hand we can compute the mag- 
netostatic field H(p, z) at any point. Also, the uniform 
magnetization state is recovered with 8 = 0, where the 



(4) functions reduces to 
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2. Vortex domain wall 

The magnetization of the vortex wall in a nanotube can 
be expressed in a similar way to the magnetization of the 
mixed state, that is M(z) = M s sin Q(z)4>+M s cos <d(z)z, 
but now 

( for < z < z w - | 

e vw (z) = \ l (z - z w ) + § for z w - | < z < z w + % 
[ 7T for z w + | < z < L 

(10) 

as presented iniS. Here Zu, is the location of the cen- 
ter of the vortex wall, and w is the wall width. Within 
the above model, the magnetostatic field can be also ex- 
pressed with equations (7) and (8), but now the functions 
f s and f v become 

fVW ^_ £ - q \z\ _ e - q \z-L\^ (u) 



where 

M s f°° 

A (P' Z ) = —J ig{<l)Ji{<lP)^{z,q)dq, (17) 

and 

S(z,q)= [ L m x (z')e-i^ z '\dz'. (18) 
Jo 

Here, m x (z) = sin<d TW (z), with TW (z) = O vw (z). 
The function s(z,q) possesses the information about the 
DW microstructure and can be straightforwardly inte- 
grated. Finally, notice that the widths of the vortex and 
transverse walls are essentially different, and should be 
evaluated within the framework presented in 19 . 
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-q\z-y\ 



dy. (12) 



B. Transverse domain wall 

The magnetization of the transverse wall in a nanotube 
can be expressed in a similar way to the vortex wall, but 
with M(z) = Af s sin9(z)x + M s cos 9(z)z>- 9 , and 9(z) 
being the same of the vortex wall. It can be noted that 
the main difference between vortex and transverse walls 
(in nanotubes) is the direction of the in-plane magnetiza- 
tion, which makes the nature of the wall properties to be 
different. The model that we use for the transverse wall^ 
leads us to understand some basic properties, like the DW 
width or the energy barrier for the reversal process. Here 
we are interested in the visualization of the magnetostatic 
field produced by a static domain wall. It can be shown 
that this field can be written as H = — VC/(p, </), z), with- 
out the usual axial symmetry of vortex-like states. The 
case of a vortex wall presents a magnetostatic potential 
independent of the angular coordinate <p, which secures 
angular symmetry. Following^, the field can be writ- 
ten as H TW = Hj w p + H™ + H T Z W i. Provided the 
magnetostatic potential of the TW can be written as 



U 



TW 
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VW 



A(p, z) cos4>, 



(13) 



the components of the dipole field created for the TW 
can then be expressed as 



H T p w {p,^z) 



HY W (p,z) 



dA(p, z) 
dp 



■ cos 
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TW 



(p,(/),z) = -A(p,z)s,m<t>, 



Hi w (p,^z) = Hr(p,z) + 



dA(p,z) 



dz 



■ cos 



0, (14) 



(15) 



0, (16) 



III. RESULTS AND DISCUSSION 

In the above sections we have obtained explicit ex- 
pressions to calculate the magnetostatic field for different 
magnetic states in nanotubes. In what follows we present 
plots of the magnetostatic field produced by the magnetic 
states relevant to nanotubes. We begin the discussion 
with the equilibrium states (U and M) and then we ana- 
lyze the case of vortex and transverse walls propagating 
along the nanostructure. 



A. Vector plots of the equilibrium states 

As mentioned earlier, the uniform magnetization state 
(#o = 0) is a special case of the mixed state (9 > 0), and 
occurs when the tube radius is approximately less than 
iL x , provided the exchange energy dominates over the 
dipolar oneM. When R > AL X the magnetization at the 
tube ends is not entirely parallel to the cylindrical z axis, 
instead it develops vortex domains at the ends in order 
to reduce the dipolar energy^. The long range char- 
acter of the dipolar interaction makes that the dipolar 
energy grows as ./V 2 , whereas the exchange energy grows 
as N, with TV the number of atomic magnetic moments in 
the tube. Therefore, the edge vortex domains increases 
in size with the radius. In the language of the above 
equations, the size of the vortex domains is given by the 
parameters 9q and d. As R increases over AL X , 0$ and d 
increases, as the reader can see from figures 5 and 6 iiJ^. 
To calculate the magnetostatic field we need to know the 
set of parameters (#0j d) which minimize the total energy 
of the mixed state. Using the expressions for the total 
energy ofi^, we obtain for nanotubes with (3 — 0.9 and 
L = 1000-La, different values of (9o,d) depending on their 
radius. These results are summarized in table ??. Note 
that for L = lOOOLz the transition to a flux-closure vor- 
tex state occurs approximately at R s=s TJLj^. 
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R/L x 9o d/L x 
4 

6 48.7 15.5 
10 75.2 41 
20 83.3 202.6 
27 83.8 414.3 

(0o,d) for magnetic nanotubes with j3 = 0.9 and L = 
lOOOL^ as a function of their radius. 




FIG. 1. (Colour online). Magnetostatic field produced by a 
tube with uniform magnetization along the symmetry axis. 
The arrows have been normalized to the local strength of the 
field indicated by the colour. Parameters are indicated in the 
text. 

The arrows in figure Q] show the normalized magneto- 
static field for a uniformly magnetized MN with R = AL X . 
The colour code indicates the strength of the field at 
each point; white colour indicates that the magneto- 
static field is stronger at the tube ends because "mag- 
netic charges" (cta/) resides only at the surfaces such 
that <jm = M • n ^ 0, with n a unit vector normal to 
the surface of the sample^. Provided that the U and M 
states (as well as the vortex wall) can be represented with 
M = MAz)<f) + M z (z)z, it follows that the surface mag- 
netic charges are given by <jm(z) — M z (z)z ■ n, where 
z = or L. In the upper surface of the tube (z = L) 
n = z and then <tm(£) = M Z (L), whereas in the lower 
surface n = — z, and ctm(0) = — M z (0). If the actual 
equilibrium state is the U or the M state, we find posi- 
tive magnetic charges (<jm(L) > 0) located at the upper 
surface, and this region behaves like a source of charges, 
as in figure [T] in elementary electromagnetism the field 
lines emanate from positive charges. Conversely, in the 
lower surface om(0) < and the region behaves like a 
sink of charge (the field lines converge to these points). 
It can be shown that for nanotubes in the mixed state, 
the dipole field looks very similar to the one plotted in 



figure [T] As a consequence, we do not present the same 
plots for tubes in the mixed state, instead we show in 
figure [5] the dependence of the magnitude of the dipole 

field (Hd = ^ Hf, + Hi ) with the radial coordinate for 

different tube radius and z = 1050-La, (just 50L X above 
the upper surface of the tube). We have normalized all 
the curves to M S (R/ L x ) 2 , and in order to remark the dif- 
ference between the U state and the M state, we present 
groups of curves for both states. 
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FIG. 2. (Colour online). Magnetostatic field strength as a 
function of the radial coordinate and z — 10507^. The left 
panel shows curves for the uniform state and the right panel 
shows the magnetostatic field for the mixed state. Different 
curves represent different tube radius, as depicted in the fig- 
ure, and parameters have been taken from table 1. 

We observe that a proper consideration of the equi- 
librium state reduces considerably the dipole field. This 
effect becomes critical as the tube radius is increased, 
which increment the size of the vortex domains of the 
mixed state**. It is worth to mention about the impor- 
tance of a proper consideration of the actual magnetic 
state in tubular nanostructures. In the sake of simplic- 
ity, it is a very common practice in magnetism to ap- 
proximate the magnetization as a uniform field, in spite 
that uniform magnetization can be achieved (at zero ap- 
plied field) only in nanostructures with dimensions of the 
order of the exchange lengt h 29 ' 30 , or when the nanopar- 
ticle is an ellipsoid with the appropriate size. Moreover, 
cylindrical magnetic nanowires and nanotubes are fre- 
quently consider as infinite structures, and only in this 
limit the approximation of a uniform dipole field becomes 
plausible, because the cylindrical nanostructures can be 
considered as ellipsoids only if they are infinitely long. 
However, actual magnetic nanostructures are not infinite 
and finite size effects can be relevant. It is well known 
that the dipole field produced by a nanoparticle with an 
ellipsoidal shape is uniform if its size is below the single 
domain limit (in absence of an external field}2£. In this 
limit one can safely write the dipole field in terms of the 
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corresponding demagnetizing factor: H f ; = — 7VM. 
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FIG. 3. (Colour online). Component of the dipole field along 
the MN axis as a function of coordinate z. The fields have 
been normalized to N Z M S , where N z is the corresponding 
demagnetizing factor. The long dashed line corresponds to 
the U state for R/L x = 4. The short dashed line corresponds 
to the U state and R/L x — 10. The thick solid line correspond 
to U state and R/L x = 20 whereas the circles correspond to 
the M state for R/L x = 20. The horizontal line illustrates the 
strength of the dipole field calculated in the usual ellipsoid-like 
approximation for uniform magnetization: Hd — ~N z M a . 

The difference between uniform and non- uniform dipo- 
lar fields is illustrated in figure [3J where we have plotted 
H z , the component of the magnetostatic field along the 
tube axis, as a function of z for different tube radius and 
magnetic states. The horizontal solid line depicts the 
strength of the (uniform) dipole field (N Z M S ) produced 
by an infinite MN with uniform magnetization, whereas 
the long dash, short dash and solid lines correspond to 
the field produced by a MN with a uniform magnetiza- 
tion state with radius R/L x = 4, 10 and 20, respectively. 
The circles correspond to the dipole field produced by a 
MN with a mixed state rather the incorrect uniform mag- 
netization state, and with R/L x — 20. The fields have 
been normalized to N Z M S , where N z = N Z (R, L, (3) is the 
corresponding demagnetizing factor for a nanotub o 12 ' 31 . 

B. Vector plots for a vortex wall 

In order to illustrate the basic characteristics of the 
magnetostatic field produced by a MN with a vortex 
domain wall, we have chosen nanotubes with /? = 0.9, 
L = lOOOL^ and R — 10L X . Following a previous work 
by some of u*^, we can calculate numerically the DW 
width (w) by minimizing the total energy including ex- 
change. It is found that w = w(f3, R), and for the values 
mentioned above we obtain a wall width w = 109L X . 
In figure 2] we show the magnetostatic field normalized 
to his strength for a vortex wall located at z w — 750L X , 
whereas figure [5] shows the same plot for the VW confined 
at the middle of the tube (z w = 500-La;). 

Both figures illustrate three regions where there are 
magnetic charge; upper and lower tube ends, and the 




FIG. 4. (Colour online). Vector plot of the magnetostatic 
field for a ferromagnetic nanotube with a vortex domain wall 
located at z w — 750L X . The colours indicate the strength 
of the magnetostatic field. Parameters are presented in the 
main text. 




FIG. 5. (Colour online). Vector plot of the magnetostatic 
field for a ferromagnetic nanotube with a vortex domain wall 
located at z w = 500L X . The colours indicate the strength of 
the magnetostatic field. Parameters are presented in the text. 



region of the DW around the wall position z w , which 
has been highlighted with red colour. There are nega- 
tive surface charges at the tubes ends, because ctm(0) = 
-M z (0) = -M s < and <r M {L) = M Z {L) = -M s < 0. 
Also, there is a positive volume charge localized in the 
domain wall region. A volume magnetic charge^ is 
given by pm = — V • M , and thus in our case pu — 
—dM z /dz. Besides, M z = M s cosQ(z), and therefore, 
Pm = M s (n/w) sin 9 > 0, just in the DW region. There- 
fore, the DW acts as a source of magnetic charge. Note 
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that we are considering here head-to-head DWs, and in 
order to analyze the case of tail-to-tail walls, we should 
replace M z by — M z and the corresponding magnetic 
charges must have the opposite sign as discussed here 
for head-to-head walls. 
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FIG. 6. (Colour online). Strength of the magnetostatic field 
for a MN with a vortex domain wall located at z w /L x = 500. 
Parameters are the same of figure [5] while the dashed region 
represent the region occupied by the magnetic matter. 

Finally, in figure |6] we illustrate the strength of the 
dipolar field as a function of the distance to the z axis, 
and for different z values. We have considered the same 
parameters as in figure [5] which means that we have mag- 
netic matter just for p/L x = 9 — 10 and z/L x = — 1000. 

C. Vector plots for a transverse wall 

To visualize the dipole field produced by a transverse 
wall, we have chosen MNs with the same parameters of 
the above section, that is (3 = 0.9, L — 1000L X and 
R = 10L X . Following a previous worki£ we can calcu- 
late the transverse DW width (w) by minimizing the 
total energy, as stated before. Again, it is found that 
w = w(/3, R), and for the values mentioned above we ob- 
tain w = 7.67L X , a rather small value as compared to 
the corresponding DW width for the VW of the above 
section. The magnetic field for the transverse wall con- 
figuration is given by equations (|14H16p . Note that the 
angular component of the field is not zero (unless 
sin</> = 0) and therefore the field does not present az- 
imuthal symmetry as the magnetostatic field for the vor- 
tex DW. It can be shown that for the geometry discussed 
here, the vector field along the plane defined by (f> = lies 
completely in that plane, where H^cj) = 0) = 0. Along 
this plane it is found that the field is very similar to the 
field for the VW (see figure [5]). This can be understood 
by noting that, mathematically, the field for the TW is 
very similar to the case of the VW. The main differences 
are the wall width and the term dependent on the angle 
4>, and the function A(p, z) defined in (|17[) . 



In figures [7] and |8] we show the magnetostatic field pro- 
duced by the magnetization of a transverse DW located 
at z w = 500 L x . Figure [7] depicts the top view of the 
field along the plane z = 600^, and figure |S] shows the 
plane z — 1200i K . Both plots represents the in-plane 
component of the field (Hj w p + Hj w <fi) normalized to 
their strengths given by the colour code. The central 
ring illustrates the top view of the nanotube. We have 
also included information about H^ w which is given by 
the length of the arrows. It can be noted that near the 
ferromagnetic tube, the z component of the field is more 
important, and thus the arrows are smaller. 




FIG. 7. (Colour online). Magnetostatic field for a ferromag- 
netic nanotube with a transverse domain wall located at the 
middle of the tube (z w = 500L X ). The plot represent the top 
view of the field, at the plane z = 600L X . The colour indi- 
cates the strength of the magnetostatic field, and along the 
contours the value of the field is constant. Parameters are 
presented in the text. 



IV. FINAL REMARKS 

In summary, we have investigated the dipole field pro- 
duced by equilibrium and non-equilibrium domain wall 
states of magnetic nanotubes. Using a continuous model 
we have obtained simple expressions to evaluate the 
dipole fields. On one hand, we can conclude that the con- 
sideration of magnetic states, as the non-uniform mixed 
state instead the idealized uniform state, reduces con- 
siderably the magnetostatic field. Thus, it is important 
for researchers to have in mind this approach when they 
compare with experimental results, in order to avoid the 
overestimation of stray fields. On the other side, mag- 
netic domain walls, which can be manipulated by exter- 
nal fields or spin-currents, introduce volumetric magnetic 
charges into the nanostructure. These charges move in 
the direction of the applied field if the DW is a head- 
to-head wall, and in the opposite direction if the wall is 



7 




a tail-to-tail one. Finally, we have concluded that the 
magnetostatic field strongly depends on the mechanisms 
of magnetization reversal. Our results are intended to 
provide guidelines to use the magnetostatic field gener- 
ated by tubular nanostructures in prospective applica- 
tions such as the generation of a magnetic trap. 
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